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I. INTRODUCTION 



The interest in the nuclear responses at intermediate energies, in particular in the quasi- 
elastic peak (QEP) region, is still widespread: indeed, the body of experimental data 
presently available in this energy regime is not easily amenable to a unified and consistent 
understanding. 

For example, the maxima of the (p,n) cross-sections [1] appear hardened, over a wide 
range of momentum transfers, by about 25 MeV with respect to the ones observed in the 
(p,p') reaction [2], an outcome still waiting for a convincing interpretation in terms of either 
many-body or reaction-mechanism effects. 

The issue is further complicated by the experiments with heavier probes: in fact, the 
hardening of the QEP position found in the ( 3 He,t) [3] and in the (d,2p) [4] reactions at 
low transferred momenta turns into a softening when the momentum transfer is increased. 
Although this behaviour might be related to the composite nature of the probes, it is notable 
that the same effect it is observed also in pion single charge-exchange reactions [5]. On the 
other hand, the predicted hardening of the isovector transverse crxq response appears well 
established according to the available data of deep inelastic (e,e') scattering [6]. 

In this paper, we shall mainly concentrate on the crq and crxq isovector channels. 
Actually, some evaluations of the spin-isospin (err) responses have already been performed 
in the framework of the bound state [7-11] or continuum [12,13] random phase approximation 
(RPA) with the conventional particle-hole (ph) interaction for the or channel. As it is well 
known, the latter includes, beside the Landau-Migdal parameter g' (or, equivalently, a short 
range interaction of some sort), the exchange of the pion and the rho meson. 

A suitable test for these calculations is represented by the above mentioned data of deep 
inelastic electron scattering, where the crxq response has been separated out in a few nuclei 
over a wide range of momentum transfers. Here the virtual photon is exploring the whole 
of the nucleus, which thus has all its constituents responding to the probe. 

When, however, one comes to consider strongly interacting probes, the ones offering 
a real hope of disentangling the elusive cr-q response, the additional problem is faced of 
appropriately dealing with the reaction mechanism, or, in other words, with the distortion 
of the impinging and outgoing hadronic waves. This is essential, since here the response of 
the nucleus is mostly confined to the surface region and one has to accurately assess how 
much of the latter is actually involved in the process. 

In this connection, the authors of Ref. [13] resorted to the Distorted Wave Born Ap- 
proximation (DWBA), whereas in Refs. [10,11] the Glauber theory [14] in the one-step 
approximation was employed. The former approach, as it is well-known, is fully quantum- 
mechanical; however, in the high energy limit and for nearly forward scattering the DWBA 
practically coincides with the eikonal approximation and thus is equivalent to the Glauber 
theory: but the latter allows one to deal with multistep processes, which are likely to give a 
substantial contribution at large momentum and energy transfers. 

This is indeed what happens and we shall explore the contribution the two-step processes 
provide, particularly in charge-exchange reactions, where their importance is expected to be 
greater. 

Turning to the many-body aspect of the problem, we calculate the or nuclear responses in 
the continuum RPA with the (g' + n + p) ph interaction, but utilizing a method different from 
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the one originally proposed by Shlomo and Bertsch [15], which has been widely exploited in 
this kind of calculations. Our approach is entirely worked out in momentum rather than in 
coordinate space and allows one to incorporate in a natural way the spreading width of the 
ph states. 

In fact, continuum RPA naturally accounts for the escape width of the particle states, 
ignored, instead, in a harmonic oscillator basis. However, the spreading width of the states 
is also quite important for a realistic description of the ar responses. The inclusion of the 
latter is achieved in Ref. [13] through the direct introduction of an optical potential in the 
Schroedinger equation for the single particle (but not for the hole) states. The disadvantage 
of such an approach lies in the ensuing lack of orthogonality between the single particle wave 
functions that needs then to be cured with a rather cumbersome procedure [16]. 

Alternatively, Smith and Wambach [17] introduced a phenomeno logical complex and 
energy-dependent self-energy that directly couples the ph to the 2p2h sector of the Hilbert 
space, a procedure that should simulate the much more computer-time consuming Second 
RPA (SRPA) [18]. We adopt here an analogous approach: however, instead of introducing 
the coupling at the level of the RPA eigenstates we do so at the level of the bare ph states and 
then proceed to the construction of the RPA solution. The two approaches differ since the 
validity of the former requires the near diagonality and state independence of the "collective" 
self-energy in the RPA basis, whereas in the latter the ph self-energy in the ph basis is 
required to be almost diagonal. Notably, the two frameworks lead to quite similar results, 
as far as the nuclear response is concerned. 

In this paper we shall test the above outlined model against deep inelastic electron scat- 
tering data (volume responses) on 40 Ca as obtained in Saclay [6]. For the surface responses 
we shall analyze the data obtained with the charge exchange (p,n) reaction at Los Alamos 
[1]. We shall consider as well the (p,p') data [2]. In both cases the proton probe has an 
energy of 795 MeV. This set of data represents a large body of what is nowadays available 
on the experimental side. We do not analyze here the (d,2p) and ( 3 He,t) data (see, however, 
Ref. [11]): the internal structure of composite probes would require an appropriate treat- 
ment, since it has a strong impact on how the probe's and target's spins are coupled [19], 
and might also qualitatively influence the distortion effects [20]. 

We organize the paper as follows: in Sec. II we present our formalism of the continuum 
RPA and the way we deal with the spreading width problem; in Sec. Ill we discuss the one 
and two-step contributions to the surface response within the Glauber theory, taking into 
proper account the cylindrical geometry of the reaction. 

In Sec. IV we show how the various ingredients of the model affect the response functions 
and then test the model against the already mentioned body of experimental data. 

Finally, in the concluding Section we summarize our results. 

II. VOLUME RESPONSE FUNCTIONS 

A. The Polarization Propagator and the Continuum RPA 

As it is well-known, the nuclear response function to an external probe is obtained 
through the imaginary part of the polarization propagator [21] (non-diagonal in momentum 
space for a finite system), which reads 
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n(q,q» = ]T < ^ |O(q)|^n >< 6^)1^0 > 

1 1 



X 



where 



huj - (E n - Eq) + irj Tiuj + (E n - E ) - irj 



H\ip n >= E n \ip n >, 



(2.1) 



(2.2) 



H being the full nuclear Hamiltonian and 0(q) the second quantized expression of the vertex 
operator. We shall mainly be concerned with the ar nuclear excitations: in such a case one 
has 



L (q,r) = T a <r-qe^ 



(2.3a) 



in the longitudinal channel and 



T (q,r) = -|<TXqe^ r 



(2.3b) 



in the transverse one. In the following, for sake of convenience, we shall set E = 0. 

The angular part of II(q, q'; uj) can be handled through a multipole decomposition which 
reads [22] 



n L (q,q';cu) = "£Uj(q,q';u)Yj; M (q)Yj M (q!) 

JM 



(2.4) 



in the crq channel and 



2 J + 1 

n T (q,q / ;w)= ^ U JJ'(Q, g'i ^jV(q)*J'M'(q') 9 v - ; 

JJ'M' ZJ + 1 



(2.5) 



in the crxq one. 

The latter expression, somewhat more involved than the former, is shortly derived in 
Appendix A. Also 



nj(g, g'; w) = pj(g, q'; 



[2.6a) 



and 



n^g,^) = ^E[nj(?,9 , ;^)] // ,6j>i> 



(2.6b) 



where 



£ 1 J 




^2.7a) 



and 
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bj' = V6V2J' + 1 a.j, e 



J J' 1 
1 1 £ 



(2.7b) 



From (2.4) and (2.5) the following expressions for the or nuclear responses are then 
obtained: 



Rl,t(q, uj) = ImIT LiT (g, q, ; u) 



with 



7T 



4tt 2 



Im^(2J + l)U J{LtT) (q,q;u), 



U J(L) (q,q';uj) = Uj(q,q';u) 

= E[nj(?,?';w: 



(2.8) 



(2.9a) 



and 



n jit) (g, q';uj) = J2 n j j' q'\ 



1 2 [Hj(g, g'; w)l (<W - a^a^)- 



(2.9b) 



The evaluation of the J-th multipole ITj, the dynamical part of the polarization propagator, 
is carried out in this paper in the RPA framework by solving the following set of coupled 
(by the tensor interaction) integral equations [22] 



n°j(g,g> 
i r°° 



(:(" 



+ 



^ / dkk 2 J2[nj(q,k;u;)} [Uj(k,u)] tlta Ilf A (M> 



(2.10) 



The bare ph polarization propagator (obtained by replacing the full nuclear hamiltonian H 
with the mean field one H°) reads 



n°,M»L „ = £<$(?) 



^ - (e p - e h ) +irj hu + (e p - e h ) - «?y 



Q^'V), (2.11a) 



where 



Q J P l(<l) = <J P Jk, J\ta; J> 5 aA (-tY +1 (-lYH[n(2£ p + l)(2£ h + l)] 1 ^ 2 



x I f 



£n p E p j p n h e h j h 



(q) r p h • 

W j 



(2.11b) 
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In Eq. (2.11b) <j p jh', J\icr; J> is the standard LS — jj recoupling coefficient and 1 

Zen p e p j p n h e hjh (q) = J o drr 2 j e (qr)R ripepjp (r)R nhehjh (r). (2.11c) 

If the spin-orbit term of the nuclear mean field is neglected, then (2.11a) is diagonal in the 
orbital angular momentum, i. e. [tl j(q,q';u)}ui = 5wll%q,q';uj), where 



fl° e (q, q'; u) = 16tt ]T (2l p + l)(2£ h + 1) 

1 



n h e h 







Ztn p i p n h i h (q) Zen p e p n h e h (q ) 



x 



huj - (e np e p 



t h )+irj hw + (e np£p - e nh e h ) - 



(2.12) 



In (2.11a), (2.11b) and (2.11c) R p (h) are the radial wave functions and e p (M the associated 
eigenvalues. They are obtained by solving the Schroedinger equation with the Woods-Saxon 
potential 



W(r) = 



1 _)_ e (r-R)/a 



he 

mlc 2 



,(r-R)/a 



ar [l + e (r-R)/ a y 



£■ a 



(2.13) 



(in the present paper, for sake of simplicity, the Coulomb potential is neglected), where 
m n is the pion mass. Obviously, for the states in the continuum part of the spectrum the 
sum over n p has to be changed into an integral, since the principal quantum number n p 
becomes a continuum variable. Instead of following the standard procedure of calculating 
the propagator in the coordinate space and then Fourier transform to the momentum space, 
we have chosen to evaluate directly Eq. (2.11a): the calculation is fast and straightforward 
and one only needs to take some care in performing the integrals over the ph energies in the 
resonance region. 

As already mentioned, in the RPA equations we employ the ph interaction g' + V n + V p , 
namely 



[Ujik^)}^ = V L (k,uj)a jei a je2 + V T {k,uj){5^ 2 - a jei a J( , 2 ), 



(2.14) 



where 



9 + 



A- 2 



V T (k,u) 



i TV 



k 2 - nl 



J O 



(2.15) 



k 2 



Ulu 2 -k 2 -y 2 



1 With respect to previous work [22], we have changed here the definition of Ol,t i n (2.3a) and 
(2.3b) by using the unit versor q instead of the vector q in the spin longitudinal and transverse 
operators. Accordingly, the quantity I is now defined without a factor q in front whereas in the 
ph interaction (2.15) appears an additional factor k 2 . Also, a factor l/y/2 has been introduced 
in the (2.3b) in order to have the same normalization for both the longitudinal and transverse 
propagators. 
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In the above, f^/Anhc = 0.08, fp/fXp = 2.18(/^/74), hp = (u;, k) and the usual monopole 
form factors 

r «w - A ?;_;rr F <««) 

have been included at the 7r(p)NN vertices with cut-offs = 1300 MeV and A p = 1700 MeV, 
respectively. Moreover, in the following we shall always use g' = 0.7. 

Since we shall also be interested in (p,p') reactions, we have to introduce, in addition 
to the spin-isospin responses, a response function Rst in the scalar (S = 0) and isoscalar 
(T = 0) channel ((p,p') scattering at intermediate energies is in fact largely dominated by 
this channel). This response reads 

Roo(q, u) = -^ImEj(2J + l)Uj m (q, q; u), (2.17) 

where at zero-order Tlj( 00 )(q,q';uj) coincides with the [Tlj(q,q';uj)]jj of Eq. (2.11a), substi- 
tuting <5 CTi i with 5 a> o in (2.11b). In RPA also Ilj( 00 ) obeys an integral equation similar to 
(2.10) and for the ph interaction we have used the G-matrix of Ref. [23] at a density of 
roughly one-half the central nuclear density. 



B. The Spreading Width 

The continuum RPA framework naturally accounts for the escape width of the particle 
states. However, also the spreading width of the latter plays an important role in the nu- 
clear many-body problem and should be reckoned with. For this purpose, it is convenient, 
following Ref. [17], to recast the expression (2.1) for the polarization propagator as follows 

II(q,q» =< ip \6(q)G(u) O f (q')|^o >, (2.18) 

where 

G(u) = 1 1 (2.19) 

%uj — H + irj hu + H — irj 

is the propagator of the excitation induced on the exact nuclear ground state \ipo > by the 
operator O(q). 

Now, in the framework of Feshbach's formalism, starting from 

H = H + V (2.20) 
one can derive the following effective Hamiltonian, restricted to operate in the ph space: 

H eff = PHP - PVQ 1 -QVP 

QH Q — Tluj — ii] 

= PHP - U(lo) 

= H° eff + PVP, (2.21) 

having set 
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H° eff = PH P-U{uj). (2.22) 

In the above P and Q are the projection operators in the ph and in the complement space, 
respectively. 

Therefore, to allow the formalism to encompass the coupling of the ph excitations to the 
2p2h ones (or to more complicated configurations), i. e. to incorporate the spreading width 
of the ph states, one should replace in the previously outlined RPA scheme the "bare" ph 
propagator with the expression 

n°(q,q';ccO = ]T < o |6(q)|<V >< ^^(q')!^ > 

ph 

1 1 



X 



hu + H p h(u) - (e p - e h ) +ir) hu - T, ph (-uj) + (e p - e h ) - ir] 



[2.23) 



where 



and 



E pft (u;) =< <j> ph \VQ(QHQ -hu- iv^QV^ph > (2.24) 



H \(j) ph >= (e p - e h )\<p ph >, (2.25) 

the \4> p h > obviously being the "bare" ph states. 

A microscopic calculation of I] ph (uj) is difficult to carry out. An alternative approach 
[17], to which we shall adhere, assumes that the real and the imaginary parts of T> ph , namely 

Z ph (uj)=A ph (u J ) + i^^, (2.26) 

can be cast into the form 

F p h(uj) = 7 P (hiu + e h ) + 7ft(e p - hu) 

(2.27) 

Ap A (cj) = A p (hu + e h ) + A h (e p - hu), 

where the arguments of the functions appearing in the right hand side of the above ex- 
pressions are inferred from the analysis of the second order diagrams of Fig. 1(a). Note 
that the diagrams corresponding to the exchange of a bubble between two fermionic lines, 
conjectured to be small in the or channel, are neglected (Fig. 1(b)). 
Now, instead of calculating these diagrams, we utilize the formulae 

(2.28) 

7h(e) = 7 P (-e), 
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symmetrical with respect to the Fermi energy (e^ = 0), which give a reasonable fit of the 
particle widths for medium-heavy nuclei, using a = 10.75 MeV, eo = 18 MeV and e\ — 110 
MeV [24]. We then insert these expressions into the once-subtracted dispersion relation 



AplCj 2tt Jo (e'-e)e'' 



(2.29) 



(and in a similar one for A h ) obtaining 
1 ae\ ( ( 



A P (e) 



7i e\ - e\ 



e 2 + e 2 e 2 + e 



— e 



log |e/e | log |e/ei| 



e 2 + e 



e 2 + e 



(2.30) 



A h (e) = -A p (-e). 



The subtraction at the Fermi surface avoids the double counting of the smooth background 
in the single particle energy already embodied in the mean field. 

Formulae (2.28) and (2.30) are the ones we shall employ in analyzing the ar nuclear 
response functions, dressing the bare ph propagator and then solving the RPA equations. 



III. SURFACE RESPONSE FUNCTIONS 
A. The Response Function in the One-Step Glauber Theory 

In order to treat the response of the nucleus to a strongly interacting probe, one should 
enlarge the framework described in the previous Section to account for both the distortion 
and the absorption of the impinging probe on the surface of the nucleus. 

In the framework of Glauber theory [14] , the scattering amplitude of a probe on a nucleus 
of mass number A is given by 

FfM) = ^J dhe^ h <^/|r(b; Sl ...s A )|^>, (3.1) 
where the so-called nuclear profile function T is expressed as 

r(b;s 1 ...s A ) = l-n[l-r j (b-s J )] (3.2) 

i=i 

in terms of the single nucleon profile function 

V ^ h) = ^nkS dXe ~ lX hf ^ X) - (3 ' 3) 

In the above, b is the impact parameter and q = k — k' the transferred momentum (k and k' 
are the projectile incoming and outgoing momenta, respectively): they are bi-dimensional 
vectors in the plane orthogonal to the direction of motion of the projectile. The probe- 
nucleon amplitudes /(A) of Eq. (3.3) are assumed to be the free ones and are meant to be 
evaluated in the laboratory system. 
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If the excitation energy of the nucleus is supplied by the probe in a single collision 
(one-step assumption), then, for large enough nuclei, one can rewrite Eq. (3.1) in the form 

FfM) = ^J db e *V**<*> <^f\T aT (b - sM>, (3.4) 

where, for definiteness, we consider the spin-isospin inelastic channel. The complex phase 
shift Xopt, responsible for the absorption and distortion of the probe, is given by 

Xo P t(6) = y/(0)T(6), (3.5) 

with 

f+OO 

T(b)= dzp{r = Vb 2 + z 2 ), (3.6) 

J — oo 

p(r) being the nuclear density (assumed in the following to be represented by a Fermi 
distribution) and /(0) the forward total NN scattering amplitude. At high energies, the 
imaginary part of /(0) is dominant and, therefore, making use of the optical theorem, one 
can write 

Xo P t(&) = ^totT(6), (3.7) 

where <7 to t is the effective probe-nucleon total cross-section (effective because empirically 
embodying Pauli blocking effects). 

Eq. (3.4) suggests to replace the vertex operators LiT (q, r) of (2.3a) and (2.3b) with 
the new ones 

°^ /(q ' r) = (27r)4r(g) / rfbrfAeiX0Pt(b)et(q_A) ' b ^^( A )°^^ A ' r )' ( 3 - 8 ) 

where the /l,t(A) are the elementary isovector spin-longitudinal and spin-transverse probe 
nucleon scattering amplitudes. From the above expression one sees that the probe does 
not transfer a single momentum q to the nucleus, but rather all possible momenta A with 
weight /l,t(A). Furthermore, its distortion is controlled by the factor exp[ix opt (6)]. The 
normalization in Eq. (3.8) has been chosen in order to recover the standard vertex operators 
in the limit of no distortion: 

l,t — ► L)T . (3.9) 

For sake of simplicity, we have not explicitly expressed the dependence of O™^ on the 
spin-isospin operators of the probe (see, however, Appendix B). If we insert the perti- 
nent surface vertex operators in Eq. (2.1) and take the appropriate matrix elements of 
the spin and isospin probe operators, we can finally define two new polarization propa- 
gators n^p(q, q';u;). The associated one-step response functions, Ri^ ur ^ (q,oo), can then 
be obtained by substituting, in Eq. (2.8), Uj( LjT )(q, q; w) with the corresponding surface 
expressions, U s J^[ T ^(q,q;uj), given by 
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n j"i)(g, q; w) = n J(L) (g, q; u) 



\h(q)\ 2 Jo 

and 

n j"?)(g, q; d) = n J(T) (g, q; u) 
4 



| roo poo 

+ JhW)T 2 Jo dXX l rfA,A ' Re [^( A )^( A/ )^(A,A';g)]n J(L) (A,A» 
1 r°° 

-2— — - ^ rfAARe[/2(g)/ L (A)/J J (A;g)]n J(L) (g,A;cu) (3.10a) 



{ /'(X) /"OO 

rrj^ / d\\ rfA / A'$:Re[/;(A)/ T (A')G j; (A,A';g)]n JJ ,(A,A» 

I roo 

-2^— / rfA A X: Re[/r(?) /t(A) #/'(A; g)]n JJf (g, A; w). (3.10b) 



j 

In Eqs. (3.10a) and (3.10b) we have set 

Gj(X, A'; q) =J2 c J(-m9* m {\ q)gm(X, q) 



(3.11) 



#j(A; g) = c Jtm9m{\ q), 

Em 



where 



roo 

g m (\q)= dbb{l-exp[i Xop t(b)}}J m (Xb)J m (qb) (3.12) 

J 

and 



2 (£-m- l)\\(£ + m- 1)!! 



(£ + m)!!(£-m)!! 



(3.13) 



0, £ + m odd 

1, £ + m even 



Iljj/ has been defined in Eq. (2.6b) or (2.9b). In Appendix B, we shortly sketch the derivation 
of Eqs. (3.10a) and (3.10b). 

Again, also in the scalar-isoscalar channel, one can introduce a surface propagator II : 
it is easily verified that its expression is identical to Eq. (3.10a), replacing everywhere (L) 
with (00) (/oo is the S = 0, T = NN amplitude) and setting a 2 je -> ^ in (3.13). 



B. The Response Function in the Two-Step Glauber Theory 

Multiple scattering contributions to the nuclear response function can be calculated along 
lines similar to the one-step case. However, the problem gets numerically rather involved 
and therefore we resort to an approximation that is often employed in Glauber calculations. 
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In this approximation, each contribution from the multiple scattering series to the re- 
sponse function is expressed as a volume response function times a distortion factor [25,26]. 
The latter is independent of the transferred energy and momentum and should embody the 
distortion effects. Thus, for instance, the one-step response can be written as 

R%{q,u)=V 1 R L , T {q,u), (3.14) 

where R LjT (q,uj) is the volume response of Eq. (2.8) and 

T>i = N$ = J dbT{b)e' a ^ b \ (3.15) 

iVgg being the effective number of nucleons participating in the single collision. 

A word of caution is in order here: expression (3.14) is able to reproduce the gross 
features of the nuclear response function (say, its size), but yields the same shape of the 
volume response. This is at variance, as we shall see in the next Section (see also Ref. [11]), 
with the full Glauber calculation of the previous Subsection, since there the response is also 
reshaped, a feature that cannot be overlooked, if one has to disentangle genuine nuclear 
correlations from distortion effects. This may be also the case, of course, when one comes 
to multistep processes. Anyway, the size of their contribution, although not negligible, is 
considerably smaller than the one arising from the one-step scattering, thus reducing the 
sensitivity of the response to the details of their shape. This assumption has, however, to 
be assessed (see next Section). 

For charge-exchange reactions at intermediate energies the two-step response function 
can then be defined as [26] 

4 2) t(^) = ^rrj-^yp / dc( [ du'\f L , T { q ')\ 2 R L , T { q ',u') 

x |/oo(|q - q1)| 2 i?oo(|q - q'U - cu'), (3.16) 
where q' is again a bi-dimensional vector, k is the momentum of the projectile and 



= ^JdbT 2 (b)e-^ T ^ (3.17) 



is connected to the effective number of pairs participating in the double scattering according 
to 

N ^ = (2)^(?tot)/© 2 (0). (3.18) 

In Eq. (3.16), the charge-exchange reaction, driven by the spin-isospin amplitudes, can occur 
only once and the second scattering is driven by the scalar-isoscalar amplitudes. The factor 
2 comes from the two possible orderings of the reaction. 

Non-charge-exchange reactions, on the other hand, are dominated by the scalar-isoscalar 
channel, leading to the following definition: 

Roo(q,v) = ^ ) f dq f du'\f 00 (q')\ 2 R 00 (q',uj') 

k l/oo(<?)r J Jo 

x |/oo(|q " q'l)| 2 ^oo(|q - q V - (3.19) 
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Note that by inspecting (3.16) and (3.19) one can already predict that two-step contributions 
will be twice more important in charge-exchange reactions than in non-charge-exchange ones. 
Indeed, the fact that in R^ T one of the two rescatterings is driven by |/oo| 2 makes i?^ T and 

(2) 

R 00 practically of the same size, apart from the factor 2 due to the two orderings of the 
charge-exchange reaction. 



IV. RESULTS 

Let us start by discussing the volume responses, relevant to electron scattering. We have 
performed the calculations for 12 C and 40 Ca, using for the Woods-Saxon potential (2.13) the 
following set of parameters: 

W = -54.8 MeV, W so = -lOMeV, 
R = 1.27 A 1 ' 2 fm, a = 0.67 fm. ^ ' 

In Fig. 2 we display the ar longitudinal and transverse response functions of 40 Ca at two 
values of the momentum transfer q, comparing the free and RPA responses with and without 
inclusion of the spreading width of the ph states. When the spreading width is included, 
one observes a sizable damping of the QEP and of the low-energy resonances, together with 
a broadening of their width (since the energy sum rule has to be conserved). The effect is 
essentially the same in the free as in the RPA responses. 

Note, as already mentioned in the Introduction, that in spite of the different approxima- 
tions for the inclusion of the collisional damping of the single-particle motion underlying our 
approach and the one of Ref. [17], the two calculations are in fact in substantial agreement. 
Note also that these correlations, related to the single-particle motion, affect the position of 
the QEP much less than the RPA ones. 

The transverse response function enters directly into the expression for the nuclear trans- 
verse structure function, measured in electron scattering experiments. Indeed, one has 

Sr(q,u) = ^ p -^) 2 G 2 M (ql)R T ( q ,u;), (4.2) 
/i being the nuclear Bohr magneton, fi p = 2.79, /i n = —1.91 and 

Gm( ^ } = [l + (q2-c,7^)/18.1fnr 2 p (43) 

the usual electromagnetic 7NN form factor. In (4.2) the small isoscalar contribution has 
been neglected. 

In Fig. 3 we compare the calculations (with inclusion of the spreading width) of the 
previous figure to the experimental data [6]. It is quite clear that the RPA correlations 
are successful in bringing the QEP position to the right place, but they miss the correct 
g-dependence of the strength. It is unlikely that this shortcoming be due to the specific 
model we employ for the ph interaction: in our treatment the Landau-Migdal parameter g' 
incorporates the effect of the exchange diagrams in the RPA series and this approximation 
is known to work well in the spin-isospin channel [27]. 
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Thus, to improve the accord with experiment, further contributions to the response 
function should likely be looked for beyond the ph Hilbert space of RPA, such as 2p-2h and 
meson exchange current terms. This is further suggested by the fact that in the high energy 
region of the response some strength is missing even at high momenta, where the strength 
of the QEP is fairly close to the data. See, for instance, Ref. [28] for a treatment of 2p-2h 
contributions in nuclear matter, where, indeed, a satisfactory description of (e,e') data is 
found. 

Let us now turn to the main topic of this work, namely to the reactions involving strongly 
interacting probes. In this case, a new feature of the response functions is related to the 
strong distortion of the projectile, whose strength is set by the effective total probe-nucleon 
cross-section of (3.7). 

In Figs. 4 and 5 one can see a comparison of the one-step longitudinal and transverse 
volume (no distortion) and surface (distorted) responses for 12 C and 40 Ca at two transferred 
momenta, q = 1.54 fm -1 and q = 2.31 fm -1 . In each plot besides the volume response 
(5tot = 0), two surface cases (<7 tot = 30 mb and a tot = 40 mb) are displayed, for both the 
free and the RPA responses. As it may be expected, one observes a strong quenching of 
the surface responses, together with a reduction of the importance of the RPA correlations, 
since the reaction is now mainly confined at the low density, peripheral region of the nucleus. 

Furthermore, another effect shows up, of great importance for the interpretation of the 
experimental data, namely a shift of the QEP position of the surface responses with respect 
to the volume ones: it appears to be sizable (due to the fact that RPA correlations are 
damped, it is much larger than the shift induced by the latter) and practically independent 
of <7 to t in any realistic range of values for this parameter. It also turns out to be independent 
of the elementary NN amplitudes (here we employ the parameterization of Bugg and Wilkin 
[29]). 

However, comparing the results for 12 C and 40 Ca, one also observes a shell dependence 
for this effect: this is better illustrated in Fig. 6, where the surface longitudinal response 
functions for 12 C and 40 Ca are directly compared for the case a to t = 40 mb at q = 1.54 
fm -1 and q = 2.31 fm -1 , both with and without inclusion of the spreading width. In order 
to bring the responses at the same scale, R L has been divided by N cS (see (3.15)), with 
7V eff ( 12 C) = 3.2 and A" cff ( 12 Ca) = 5.6 , respectively. 

In Ref. [11] it has been shown that the 40 Ca one-step surface responses are hardened for 
q ^$ 1.8 fm -1 and softened for q ^ 1.8 fm -1 ; in contrast, 12 C surface responses are always 
hardened, as one can see from Fig. 6. Furthermore, this feature is independent of the specific 
value of £r tot , which is dictated by the energy and the nature of the projectile. 

A contribution that might, in principle, reshape the quasi-elastic response functions is 
the two-step term in the Glauber multiple scattering expansion. In Ref. [26] it had been 
shown to produce at small transferred momenta a contribution smoothly increasing with 
the transferred energy. However, this feature, because of the NN amplitudes entering into 
(3.16) and (3.19), is strongly dependent upon the momentum regime. 

To figure it out, we plot in Fig. 7 the two-step response for a simple model, namely 
assuming the two squared amplitudes in (3.19) to be equal and with a gaussian shape, 
\f{ ( l)\ 2 = ^4exp {—f]q 2 )- Here, and also in the following calculations with realistic NN am- 
plitudes, in order to estimate we use for simplicity the free harmonic oscillator model 
without inclusion of the spreading width: since R^ is the convolution of two ph response 



14 



functions, it should be rather independent of the details of their shape. 

The curves in Fig. 7 correspond to rj ranging from 0.001 fm 2 to 1 fm 2 and it is quite 
apparent that the shape changes when the g-dependence of the amplitudes becomes more 
pronounced. The reason for this behaviour lies in (3.19): indeed, when f(q') is a rapidly 
decreasing function, the main contribution to the integrals comes from the region q' ~ and, 
as a consequence, u/ ~ 0; it then follows that at fixed q the maximum of R {2) (q,uj) will be 
found for uj around the QEP position. Note that for 77 ^ 0.1 fm 2 one gets, at high momenta, 
a maximum below the QEP and that realistic values for 77, which fit the NN amplitudes, fall 
just in this range. 

Let us now see how our full model (namely one-step RPA Glauber responses with spread- 
ing width plus two-step free responses) compares to the available (p,n) data [1]. The double 
differential cross-section for a (p,n) reaction is given by 

d2<T E \f a (q)\ 2 [R^(q,u,)+R£\qM} (4-4) 



(Kldu , , 

a=L,l 

and it is tested in Fig. 8 on the data from a reaction on 12 C at a proton energy of 795 MeV 
and for four different scattering angles (0 = 9°, 12°, 15° and 18°), corresponding to q ranging 
from 1.16 fm" 1 to 2.31 fm" 1 . 

The dashed and the solid curves in each plot represent the one-step and the full calcu- 
lations, respectively, whereas the dotted curve is the two-step term; for sake of comparison, 
we display also the cross-section corresponding to the simpler model based on the one-step 
contribution of Eq. (3.14) (dot-dashed line), i. e. 



d 2 



a 



E \f a (q)\W(q,u)+RM(q,u;)}. (4.5) 



dVLdoo , ,, 

a=L,l 

A few comments are in order. First, the reduction in strength in going from (4.4) to (4.5) 
is due to RPA correlations (indeed, the difference fades away at large q), whereas the shift 
in the QEP position comes from the cylindrical geometry of the reaction, which is respected 
by our approach (it stays constant with q). Looking at the data, it is apparent that now 
the situation is complementary to what we have observed in the case of the (e,e') data: in 
fact, we are now able to reproduce the correct g-dependence of the strength, but not the 
QEP position of the (p,n) data, which is still somewhat more hardened than predicted by 
our calculations. 

Note that the height of the peak depends on er tot : here, we have used <7 tot = 40 mb, 
without accounting for medium effects on this parameter. The estimate of this effective 
parameter is not well assessed, since there are large discrepancies between the direct cal- 
culation of Pauli blocking effects [25] and the derivation of <7 to t from the optical potential 
[26]. Using <7 to t = 30 mb, the curves in Fig. 8 would get a 25% increasing at any momentum 
transfer, without affecting the g-dependence of the strength and reproducing correctly also 
the height of the peak. 

The real problem with the model we are discussing is related to the A-dependence of 
the QEP position that it introduces (see Fig. 6), since the (p,n) data of Ref. [1] appear 
to scale, at least in 12 C and Pb. A scale factor is also the only difference in the ( 3 He,t) 
cross-sections on 12 C and 40 Ca [3], whose QEP position displays a pattern, as a function of 
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q, quite different from the (p,n) case, in agreement with our calculations for 40 Ca [11], but 
not for 12 C. In this connection, we remind the reader that in Ref. [11] the 3 He projectile 
has been treated as structureless, whereas, as already mentioned, a proper account of its 
internal structure might affect the shape of the response functions [20]. 

It is not clear, at the moment, whether this shell dependence of the Glauber quasi- 
elastic responses is a genuine effect or, rather, it reflects a shortcoming of the model or the 
approximations of our treatment: indeed, in the two-step contribution we have neglected 
the cylindrical geometry of the reaction, which, as we noted, is the source of the shift in the 
one-step term. Actually, the discrepancy between 12 C and 40 Ca shows up at large momenta 
(see Fig. 6), where the two-step response is relatively more important. 

Finally, in Fig. 9 we show the same comparison as in Fig. 8, but for (p,p') scattering at 
795 MeV [2]. We have considered only the dominant scalar-isoscalar channel, which means 
that in (4.4) and (4.5) we have set a = 00. By inspecting Fig. 9 a few observations follow: 
the two-step contribution is smaller than in the (p,n) case, as anticipated in Sec. Ill B; the 
simple calculation based on (3.14) of course yields a better account of the QEP position, 
since (p,p') data do not show any shift; on the other hand, the g-dependence of the height 
of the peak is unsatisfactorily described, but this is due to our choice of only allowing the 
scalar-isoscalar channel: when the momentum transfer increases (say, over 1.5 fm -1 ) other 
components of the NN amplitudes become important and should help in increasing the 
height of the response. 

V. CONCLUSIONS 

The aim of this paper has been to investigate the quasi-elastic nuclear response as seen in 
electron and, especially, proton inclusive scattering. We restricted ourselves to unpolarized 
scattering, since already in this domain a number of issues wait to be clarified. 

Our attention has been mainly directed to the (p,n) and (p,p') reactions, owing to the 
surprising and apparently contradictory figures displayed by these experiments. We have 
first applied our formalism, based on continuum RPA plus spreading width of the ph states, 
to the transverse electron scattering finding that: 

i) RPA correlations shift the QEP position to higher energy; 

ii) the g-dependence of the strength in the QEP domain is not correctly reproduced by 
our model, since some strength is actually missing at low but not at high momenta; 

iii) in the high energy tail of the response strength is missing at all momenta. 

The solution of the problem in items ii) and iii) likely comes from the contribution from 
2p-2h processes and meson-exchange currents, as can be seen from the results of Ref. [28], 
worked out in nuclear matter. Neglecting 2p-2h contributions in hadron scattering should 
be a less serious shortcoming, because of the strong density dependence, which makes them 
relatively less important in the low density surface regions probed by strongly interacting 
projectiles. 

In hadron scattering, on the other hand, multiple-scattering contributions can be sub- 
stantial and, accordingly, we have adopted the framework of Glauber theory: one-step pro- 
cesses have been calculated consistently within the theory, whereas for two-step processes 
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we have resorted to the approximation of accounting for the effect of distortion through a 
multiplicative factor. From the analysis of the hadron responses we have established the 
following: 

i) RPA correlations are strongly damped, but not completely washed out: this is at vari- 
ance with the polarized (p,rt) experiment in Los Alamos at q ~ 1.72 fm _1 [30], which 
seems to be compatible with free responses. It remains to be understood which mecha- 
nism is responsible of the complete suppression of the RPA correlations experimentally 
found; 

ii) in contrast to (e,e') scattering, the g-dependence of the strength is well reproduced (of 
course, properly accounting for the non-scalar-isoscalar channels in the (p,p') reaction 
at high momenta), confirming the validity of the assumption that 2p-2h contributions 
should be negligible for hadronic reactions; 

iii) the height of the QEP depends on the value of cf tot , which is not well assessed: note 
that if one uses the free value, 40 mb, the height of the (p,p') cross-section is well 
reproduced, whereas the height of the (p,n) one is underestimated; 

iv) again in contrast to (e,e') scattering, the QEP comes out in the wrong position: the 
calculations in 12 C predict a hardening of the cross-section in the whole range of 
momenta explored (~ 1 4- 3 fm _1 ), as it happens in the (p,n) but not in the (p,p') 
reaction. Items iii) and iv) together apparently point to the need for contributions 
that appear only in (p,n) and not in (p,p') scattering; 

v) the shift of the peak in the calculated responses is a consequence of the distortion of 
the proton wave and it is not the same in 12 C and 40 Ca, whereas the available data 
appear to scale. Since the Glauber theory is equivalent to the distorted wave Born 
approximation in the eikonal limit (which might or might not be valid at the incident 
energy of 800 MeV), it would be interesting to see if also in the DWBA framework a 
shell dependence shows up. Note also that the two-step response is not treated in a 
way fully consistent with the Glauber theory, a fact that could have some influence on 
the shape of the response functions; 

vi) at high transferred momenta the two-step term is sizable in the (p,n) reaction and, 
because of the momentum dependence of the NN amplitudes, does not give a smoothly 
increasing background, as it does at small momenta, but rather a contribution peaked 
slightly below the QEP. 

No final statement can be made, at the moment, on these issues: however, we wish to stress 
the necessity of formulating a theory of the nuclear response able to cope simultaneously with 
all the different kinds of reactions. Accurate calculations are surely needed, but only from 
a careful cross-referencing of the many phenomena unveiled in the scattering on complex 
nuclei can one hope to obtain a guide towards a solution of the difficult many-body nuclear 
problem. 
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APPENDIX A: THE TRANSVERSE POLARIZATION PROPAGATOR 

The transverse polarization propagator is defined as the trace of the current-current 
propagator, namely 

II T (q, q'; u) = E n„ n (q, q'; u), (Al) 

n 

where H mn is obtained by substituting (crxq) m and (crxq)„ for the two vertex operators in 
Eq. (2.1) [22]; m and n are spherical indices. 

Introducing vector spherical harmonics, H m n can be recast in the following form: 

II mn (q, q'; u>) = E n J;JlJ2 (g, q'; ^j^(q) y£g M (q'), (A2) 

JM 
•h J 2 

where 

Uj, JlJ2 ( q , q ';oo) = W [tlj(q,q';u;)] b J j\b%. (A3) 

The quantities in (A3) have been defined in Sec. II A. Using the addition theorem for the 
vector spherical harmonics, namely 

^ E E Y jTm^) ySmW) = Sj lJ2 (2J + l^q • q'), (A4) 



M=—J m 



one finds 



n T (q, q'; W ) = E ^^W^ & ^ ^ Pj '^ ' ^ 
JJ' 47r 

2 7-1-1 

= E n^^^w^^M^qOoiTXT, (A5) 

i. e., Eq. (2.5) with IT jj,(g, g'; u;) = Uj.jfj>(q, q'; u). 

Using Eq. (A3) the diagonal part of the transverse propagator can be written as 

n T (q, q;^) = E ^ E N?, E K (A6) 
j 47r «' j' 

It is easy to check that J2j'bj'ebj' e , = See' — ajeo-JV (aje being defined in Eq. (2.7a)), from 
which the result of Eq. (2.9b) follows immediately. 
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APPENDIX B: THE SURFACE POLARIZATION PROPAGATOR 



Here we briefly sketch the derivation of the longitudinal surface polarization propagator 
(3.10a). Similar calculations apply for the transverse one. 

As mentioned in Sec. Ill A, the expression (3.8) for the surface vertex operator 0™ r ^(q, r) 
is oversimplified: actually, in (3.8) one should add, inside the momentum integral, the matrix 
element of the probe spin longitudinal operator, < Sj|cr( p )-A|sj > (the treatment of the 
isospin is trivial). Inserting the resulting expression in the formula (2.1) for the polarization 
propagator, averaging over the initial probe spin and summing over the final spin Sf one 
finds (for simplicity, we consider only the diagonal part of IT L , i. e. q=q') 

nr / (q,q;o;) = * 1 / dbdX J db' dX' X ■ X' /2(A)e~ ix °p t(6) e~ i(q_A) ' b 



xIT L (A, A'; uj)f L {X')e l ^ {b ' ) e t ^' x '> h ' 



(**)*\f L (q)\ 
n L (A,A';cc 
1 1 4tt 

(27r)4|/ L ( g )pT JM n 
xUj(\, A'; ^)/ L (A0e^( 6 \ l ^- A ')- b Y; M (A)Y i ;(A)Y JM (A0yin(A / 



Y,Y.J dhdX J db'dX'fi(X) 



e *x; pt (6) e -»(q-A)-b 



where we have used the angular momentum expansion (2.4) and set X n = ^J~AnJ%Y ln {X) . 
After some algebra, one can recast (Bl) in the following form 



nsurf I \ 
l (q,q;w) 



2J+ 1 



l/^)l 2 Vtr 2 ^ + 1 



EE 



roo roc poo roo 

x / dbbe~ lx °^ b) / db'b'e iXopt{b) / dXXf* L (X) / dX' X' f L (X')Uj(X, A» 
Jo Jo Jo Jo 

xJ7 m (q;b,\)r em (q;V,\'), (B2) 
aji being defined in Eq. (2.7a) and 



Tim{q\ b, A) = ^ fj d(f> b fj d0 A e^- A )' b YUA). 
A is a versor in the (x, y) plane, hence 



(B3) 



,7T 



Y£ m (\) = Yt m {—, 4>\) 



= (-1) T 



(2£ + l)(£-m)! 
4vr(£ + m)! 



if(0)e^ 



'2£ + V 



An 



(£-m- l)\\(e + m- 1)!! 



(£ + m)!!(£-m)H 



x 



(_l)<+r 



;-i)«+"»)/ 2 + (-1)- 



{i+m)/2 



(B4) 
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where K^ m is defined as the the quantity in brackets in (B4). 

Thus, employing a standard integral representation of the Bessel function, one finds 

F em (q; b, A) = (~l) m K em J M (qb) J M (Xb) (B5) 

and 

nr / (q,q;c) = ^E(2^+i) 

roo roo 

x dXXf* L (X) dX' X' f L (X')Uj(X, A'; lo) 
Jo Jo 

xEc* dbbe-^^J m {qb)J m {Xb) / db'b'e^ b ^J m (qb')JM'b'), (B6) 



with c J£m given by Eq. (3.13). 

Using the well-known orthogonality relation 



f°° dbbJ m {qb)J m {Xb) = (B7) 
jo A 



one can write 

8(q - A) 



(°° dbb^^J m {qb)J m {Xb) = [ \ } -g m (X,q), (B8) 
Jo A 

having defined g m (X,q) in Eq. (3.12). Substituting (B8) in (B6), it is then straightforward 
to obtain Eq. (3.10a). 
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FIGURES 



FIG. 1. ph self-energy diagrams: (a) single particle diagrams; (b) interference diagrams. 

FIG. 2. ot longitudinal and transverse response functions of 40 Ca at q = 410 MeV/c and 
q = 330 MeV/c. Free (dashed line) and RPA (solid line) responses without spreading width and 
free (dotted line) and RPA (dot-dashed line) responses with spreading width are displayed. 

FIG. 3. Structure function of 40 Ca at q = 410 MeV/c and q = 330 MeV/c. Free (dashed line) 
and RPA (solid line) contributions include the spreading width. Data are from ref. [6]. 

FIG. 4. ot longitudinal response functions of 12 C and 40 Ca at q = 1.54 mi" 1 and q = 2.31 
fm _1 . Free (dashed line) and RPA (solid line) volume responses are reported, together with free 
(dotted lines) and RPA (dot-dashed lines) surface responses for two- values of otot- <?tot = 30 mb 
(higher curves) and <r tot = 40 mb (lower curves). The spreading width is always included. The 
40 Ca surface responses have been multiplied by a factor 2. 

FIG. 5. As in Fig. 4, but for ot transverse response functions. 

FIG. 6. ot surface longitudinal response functions (divided by N e g) of 12 C and 40 Ca at q = 1.54 
fm _1 and q = 2.31 fm _1 with (SW) and without (no SW) inclusion of spreading width. In each plot 
the 12 C free (dashed line) and RPA (solid line) responses are compared to the 40 Ca free (dotted 
line) and RPA (dot-dashed line) responses. 

FIG. 7. Two-step responses of 40 Ca with gaussian amplitudes |/(<?)| 2 = Aexp {—rjq 2 ) at q = 1.4 
fm _1 and q = 2.4 fm^ 1 for various values of 77: (1) 77 = 0.001, (2) r\ = 0.01, (3) 77 = 0.1, (4) rj = 0.5, 
(5) 77 = 1; the arrows show the position of the QEP. The scale is arbitrary. 

FIG. 8. Inelastic (p,n) cross-sections on 12 C at q ranging from 1.16 fm _1 to 2.31 fm _1 . The cal- 
culation based on Eq. (4.4) is reported (solid line), together with the separate one-step (dashed line) 
and two-step (dotted line) contributions and with the calculation based on Eq. (4.5) (dot-dashed 
line). Data are from the 795 MeV experiment of ref. [1]. 

FIG. 9. As in Fig. 8, but for inelastic (p,p') cross-sections on 12 C at q ranging from 1.42 fm _1 
to 1.93 fm -1 . Data are from the 795 MeV experiment of ref. [2]. 
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